Nuclear masses: evidence of order— chaos coexistence 
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Abstract 

Shell corrections are important in the determination of nuclear ground- 
state masses and shapes. Although general arguments favor a regular single- 
particle dynamics, symmetry-breaking and the presence of chaotic layers can- 
not be excluded. The latter provide a natural framework that explains the 
observed differences between experimental and computed masses. 
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Different approaches have been developed to reproduce the systematics of the observed 
nuclear masses. Some of them are of microscopic origin, while others, more phenomeno- 
logical, are inspired from liquid drop models or Thomas-Fermi approximations. Often, the 
total energy is expressed as the sum of two terms 



where Z and are the proton and neutron numbers, respectively, and x represents a set 
of parameters that define the shape of the atomic nucleus. The first term V describes 
the bulk (macroscopic) properties of the nucleus, and contains all the contributions that 
vary smoothly with proton and neutron numbers. Its typical value is close to 8 MeV/^ 
throughout the periodic table, where A = + Z is the mass number. The second term 
f/ describes shell effects, related to the shape-dependent microscopic fluctuations of the 
nuclear wavefunctions IQ . A method to incorporate these effects in a consistent manner was 
originally proposed by Strutinsky 0. 

Global nuclear mass calculations have been pursued over the years with increasing pre- 
cision Despite the numerous parameters contained in the models, the accuracy of the 
results obtained and the quality of the predictions are impressive. Different models and po- 
tentials yield similar results, and give an accuracy of 5 x 10~^ for a medium-heavy nucleus 
whose total (binding) energy is of the order of 1000 MeV. 

One may ask whether the difference between measured and computed masses has a par- 
ticular significance, since it seems to be quite model independent. Among other possibilities, 
a natural explanation would come from many-body effects not included in the mean-field 
scheme. Our purpose here is to show that there is a very natural and appealing dynamical 
explanation. We will argue that U may be splitted into two parts, U = Ureg + Uch, and 
that the present calculated masses are able to correctly reproduce only Ureg- The two con- 
tributions Ureg and Uch originate in regular and chaotic components of the motion of the 
nucleons, respectively. Although our calculations are done within a mean-field approxima- 
tion, the final result for the fluctuations produced by the chaotic part of the motion are in 
fact of a much more general validity, and may be interpreted as arising from the residual 
interactions. 

With essentially no free parameters, this dynamical symmetry-breaking mechanism pro- 
vides a quantitative description of the basic observations on nuclear masses. First, the 
amplitude of Uch is much smaller than that of Ureg- Second, the typical size of the fluctu- 
ations Ureg havc a constant amplitude as a function of the number of nucleons. The RMS 
of Ureg IS fouud to be of about 3 MeV, which is in good agreement with the typical size of 
the difference 6U = Uexp — U between the experimental and the computed bulk properties 
of the nuclei [Q. Third, the RMS of Uch is of order 0.5 MeV, in good agreement with the 
typical size of the difference SU = Uexp — Ucaic between the experimental and the calculated 
masses including shell effects The mass-number dependence of this difference is also 

well reproduced. 

From a semiclassical point of view, shell effects are interpreted as modulations in the 
single-particle spectrum produced by the periodic orbits of the corresponding classical dy- 
namics Given the single-particle energy levels Ej{x) computed from some Hamiltonian 
H{x), the level density p{E,x) = J2j S[E — Ej{x)] is approximated by 



U{Z, N, x) = U{Z, N, x) + U{Z, N, x) , 



(1) 



p{E,x) =p{E,x)+p{E,x) . 



(2) 
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The quantity p is the average density of states, whereas the oscillating part is expressed as 

oo 

P = '^^^Ap^r{E,x)cos[rSp{E,x)/h + Vp^r] ■ (3) 

P r=l 

The sum is over all the primitive periodic orbits p (and their repetitions r) of the single- 
particle Hamiltonian. Each orbit is characterized by its action Sp, stability amplitude Ap^^, 
and Maslov index Vp.r- 

The shell correction to the nuclear mass is computed by inserting the oscilla- 
tory part of the density of states into the expression of the energy, U (x, A, T) = 
/ dE E p{E,x) f{E,iJ,,T), with / the Fermi function. In a semiclassical expansion, the 
leading order of the integral comes from the energy dependence of the action, and the de- 
pendence of the prefactors can be ignored. Setting moreover the temperature T to zero, one 
obtains 

f/(x, A) = 2n'J2Y.^ cosirSp/h + z.^,,) , (4) 

p r=l ^ ^p 

where Tp = dSp/dE are the periods of the periodic orbits. The classical functions entering 
this expression are evaluated at the Fermi energy Ep, related to the mass number and shape 
parameters through the condition Jq''p{E,x) dE = A. According to Eq. (|), each periodic 
orbit produces a modulation or bunching of the single-particle states on an energy scale h/Tp. 
The presence of fluctuations in the total energy is therefore a very general phenomenon 
that occurs for an arbitrary Hamiltonian, irrespective of the nature of the corresponding 
classical dynamics. However, their importance (i.e., their amplitude) depends strongly on 
the properties of the dynamics, and in particular on its chaotic or regular character In 
order to see this explicitly, we compute the variance of the fluctuations of the energy. The 
average is performed over a mass number window around a given nucleid, as when analyzing 
experimental data. The size of the window is taken to be small with respect to macroscopic 
quantities like the Fermi energy, but large compared to the typical oscillatory scale of U. 
This average is denoted by brackets. From Eq. (^) we obtain for the variance a double 
sum over periodic orbits labeled by the indices p and p'. By considering the dependence of 
the density of periodic orbits with their period, it has been shown |]^,|^ that the dominant 
contribution to the double sum is given by the diagonal terms p = p', r = r' . The resulting 
sum is convergent. Using moreover the semiclassical definition of the form factor, 

^oo 

K{t) = 2h de cosier /h) {p{E - e/2)p{E + e/2)) 
Jo 

whose diagonal part is 

Kn{T) = h'Y.K,rH^-rrp) , (5) 

p,r 

the energy variance may be written as 
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This is the basic equation that provides the starting point for our analysis. An accurate 
evaluation of {U"^) requires the knowledge of the periodic orbits, and therefore of the Hamil- 
tonian. However, we do not want to make a calculation for a specific system. Instead we 
are interested in making general statements based on the statistical behavior of the single- 
particle orbits. That behavior depends on the regular or chaotic nature of the dynamics. 

General considerations suggest that the single-particle motion in the nucleus should be 
dominated by regular orbits. This is based on the following arguments. There are several 
factors that are important when considering the amplitude of the shell effects as given by 
Eq. (^). The first one is the prefactor Ap^r/^p- It implies that the more stable and short the 
periodic orbit is, the larger its contribution. The second is associated with the interference 
properties. Shells are wave effects: they are determined by a superposition of oscillatory 
contributions from different periodic orbits. The third factor, closely related to the previous 
one, is the phase space structure of the periodic orbits, i.e., their degeneracy. It depends 
on the existence of conserved quantities. In fully chaotic dynamics (where only the energy 
is conserved), the periodic orbits are isolated and unstable. In contrast, if there are as 
many conserved quantities as degrees of freedom, the motion is regular (or integrable) and 
periodic orbits come in degenerate families, all orbits having the same properties (action, 
stability, etc) within a family. Since many orbits contribute coherently, this degeneracy leads 
to an enhancement of the size of the shell corrections. In its equilibrium state, and for a 
given number of nucleons, the nucleus will adapt its shape in order to minimize its energy. 
This minimization procedure, that takes also into account the contributions from U, favors 
regular single-particle dynamics that produce strong effects. 

Equation (|^) allows to estimate the typical fluctuations arising from a regular dynamics. 
In the range <^ t <^ Tj^, where is the period of the shortest periodic orbit in the 
semiclassical sum (^ and = h'p is the Heisenberg time, the form factor of the regular 
single-particle levels is given by Kd{t) = t^^ |@|. This approximation is not valid for short 
times, of the order of Tmin- The reason is that for r ~ Tmin, Eq- (|§) produces system specific 
delta function peaks. If for simplicity we ignore this and extrapolate the statistical behavior 
down to r Tmin, and moreover use Kd^t) = for r < Tmin (cf. Eq. (H)), then Eq. (H) 
leads to 

The quantity Ec is the energy associated with the shortest periodic orbit, 

Ec = h/Tmin , (8) 

whereas the parameter g measures this energy in units of the single-particle mean level 
spacing 5 = 1/p 

g = Ec/S = TH/Tmin ] (9) 

g counts the number of single-particle states on the scale E^.. To compute E^ we need the 
period of the shortest orbit. Its length will typically be two or three characteristic nuclear 
dimensions. For simplicity we assume it to be two. Then, for a fiat mean-field, like the 
Woods-Saxon potential, we have Ec = irEp/kpr, where r is the nuclear radius and kp is the 
Fermi wave vector. Since r ^ I.IA^/'^ fm, we arrive at 
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(10) 



Ec is the largest oscillatory energy scale where coherent bunching effects occur in any ther- 
modynamic quantity. This result should be contrasted to the traditional shell-effect estimate 
based on the harmonic oscillator potential, Ec = hwo ~ AO/A^^^ MeV 
Since 6 ^ 2Ep/3A = 25 /A MeV, Eq. (D leads to 

9 = vrA^/l (11) 

From these results we obtain for the typical fluctuations areg = \J {U'^eg) ^ nucleus with a 
regular single-particle dynamics the value 

areg = 2.84 McV , (12) 

an expression independent of the number of nucleons. This estimate is in good agreement 
with the ~ 3 MeV observed when U is subtracted out from the experimental or calculated 
values. The use of the real nuclear mean level spacing 6 to compute ([TT|) (i.e., including 
spin and isospin degrees of freedom) amounts to treating the different contributions as 
uncorrelated, and therefore provides a lower bound. 

In deriving Eq. (|r^) a regular single-particle motion of the nucleons has been assumed. 
Even though we have stressed that a regular motion favors the minimization of the total 
energy, in the full many-body problem nothing guarantees the perfect integrability of the 
semiclassical dynamics. If some dynamical symmetries are broken, chaotic components will 
coexist with regular single-particle motion. The consequences of their presence on the 
behavior of the nuclear masses will now be discussed. 

From a semiclassical point of view, the most simple approximation that can be made 
in the generic case of a mixed dynamics, when regular orbits coexist with chaotic layers, is 
to split the sum over periodic orbits in Eq. (^) into two terms, one from the regular part, 
another from the chaotic part. The shell energy is now written as 

U = Ureg + U^H. (13) 

The two terms are, from a statistical point of view, independent, (UregUch) = 0. This 
happens because, as already pointed out, the dominant contribution to the energy comes 
from the short orbits, with Tp ^ Th- Since the orbits contributing to each term are different, 
the cross products vanish by the averaging procedure (assuming that the actions of the orbits 
are incommensurable). The decomposition (|13D identifies the difference between measured 
and computed masses, that have a typical size of ~ 0.5 MeV, with the oscillatory contribution 
of the periodic orbits lying in the chaotic components. 

We are therefore led to evaluate the variance of the shell corrections that originate in the 
chaotic layers. This can be done from the general Eq. (^). Again, our purpose is to make a 
statistical estimate valid for a generic chaotic component, with no reference to a particular 
system. In the range r^m -C t <^ r^^ the form factor of chaotic single-particle levels is given 
by Ko{t) = 2t 0. It coincides with the random matrix prediction [|l^. This result is not 
valid for times of the order of Tmin- Extending nevertheless this behavior down to r ^ Tmin, 
and imposing again Ko{t) = for r < Tmin-, we now arrive at 
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(14) 



Using Eq. (p!OD, the typical size of the fluctuations = y {U^h) of ^ chaotic component in 
the single-particle dynamics is given by 

^^'^ = ^3 MeV. (15) 

This expression has to be compared with the RMS of the difference SU between the experi- 
mental and the computed masses, given in Ref. The comparison is shown in Fig. 1. 

The agreement is extremely good. The amplitude in Eq. (|1^) is uncertain up to an 
overall factor of say, 2. It can be varied by increasing slightly the period of the shortest 
orbit (we have chosen the shortest possible one; any modification will diminish ach) and by 
the inclusion of spin and isospin (this increases ach by a factor 2 if these components are 
treated as uncorrelated) . The A dependence is very well fitted in the region A ^ 75, with 
deviations observed for lower mass numbers. This is in agreement with the limited accuracy 
of the Strutinsky corrections and in general of semiclassical theories for light nuclei. 

There are several features of the present theory that make its predictions reliable. First 
of all, Eq.(l^) contains only one physical parameter, the period of the shortest chaotic 
periodic orbit, which is a function of A because the size of the nucleus increases with the 
mass number. But it has no dependence on the relative size of the chaotic region, i.e., on 
the fraction of phase space occupied by chaotic motion. Without this quite remarkable and 
important fact we would have been forced to estimate that fraction, something that is hardly 
possible with present knowledge despite the efforts in this direction since the pioneering work 



of Ref. |TT|. This is in contrast to what happens for the regular regions. The enhancement 
factor g present in Eq. (|^) - but not in Eq. (|T^) - depends on the corresponding regular 
average level density 'p^eg- The evaluation of g in Eq. (|lT]) assumes that the nuclear single- 
particle levels are regular, and the good qualitative agreement with the experimental results 
of [/^gg suggests that indeed most of the phase space is occupied by regular trajectories. 
Mixed systems present therefore a peculiar structure of shell effects. The amplitude of the 
fluctuations of the regular phase space regions is proportional, through the factor g, to the 



square root of the regular average level density, areg oc ^Preg- contrast, the amplitude of 
the fluctuations associated with the chaotic motion is constant, independent of their phase 
space volume (aside from the transient nearly-integrable regime, not considered here, where 
the chaotic layers are formed). In the extreme case of a fully chaotic dynamics, p^g^ = and 

U = Uch- 

Equation ([l^) is in fact quite robust. It is not only independent of the chaotic phase- 
space volume, but is also valid for arbitrary dimensions. This fact suggests a very natural 
origin of these fluctuations, that we haven't discussed yet. Although our analysis is based 
on a single-particle picture, it can be extended to the full many-body phase space. In 
that space, and in a semiclassical picture, to flrst approximation the point representing the 
system follows a very simple trajectory driven by the regular mean-fleld that dominates the 
motion. On top of that, it is likely that the residual interactions, not taken into account 
in that approximation, induce chaotic motion. The presence of chaotic orbits would then 
introduce additional long-range modulations in the regular single-particle density of states. 
Eq. (|T^ would then still be valid to evaluate the amplitude of those modulations, with 
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Tmin the period of the shortest chaotic orbit in the multidimensional space. Rough estimates 
indicate that r^m is comparable to the three-dimensional non-interacting period. Therefore 
Eq. (|T5|) presumably gives an estimate of the mass fluctuations arising from neglected many- 
body effects. 

To summarize, we have investigated shell effects of nuclear ground states within a unified 
framework, namely periodic orbit theory. We show that a dynamical symmetry-breaking 
mechanism that introduces chaotic layers in the motion of the nucleons produces additional 
shell corrections to the total energy. The intensity of the effect is, to first approximation, 
independent of the size or the topology of the layers. It is governed by a single parameter, 
Ec, which is proportional to the inverse time of flight of a nucleon across the nucleus at 
Fermi energy. The comparison of the typical size of the chaotic fluctuations is in very good 
agreement with the deviations between computed and experimental values, for the amplitude 
as well as for their dependence with the mass number. In contrast, the fluctuations produced 
by the regular components are independent of the mass number, and are proportional to 
the regular phase-space volume. The picture we are suggesting is coexistence of order 
and chaos as produced, for instance, by residual interactions. It manifests in nuclear low- 
energy properties. Further evidence should be given. In particular, our theory predicts 
autocorrelations in the total energies 0], as well as the effect of the presence of chaotic 
layers on the total level density. This deserves further investigation. 

We are indebted to W. J. Swiatecki whose questions and persistence led to the present 
investigation, and to G. Bertsch and N. Pavloff for enlightening discussions and comments. 
The Laboratoire de Physique Theorique et Modeles Statistiques is an Unite de recherche de 
I'Universite Paris XI associee au CNRS. 
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FIG. 1. RMS of the difference 5U between computed and observed masses as a function of 
mass number A. Dots taken from Fig. 7 of Ref. [3]. Solid curve from Eq. (17). 
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